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Abstract 



2 Most conspicuous organisms Eire multicellulEir and most multicellulEir organ- 

3 isms develop somatic cells to perform specific, non-reproductive tasks. The 

4 ubiquity of this division of labor suggests that it is highly advantageous. In this 

5 paper, I present a model to study the evolution of specialized cells. The model 

6 allows for unicellular and multicellular organisms that may contain somatic (ter- 

7 minally differentiated) cells. Cells contribute additively to a quantitative trait. 

8 The fitness of the organism depends on this quantitative trait (via a benefit 

9 function), the size of the organism, and the number of somatic cells. This model 

10 allows one to determine when somatic cells are advantageous and to calculate 

11 the optimum number (or fraction) of reproductive cells. I show that the fraction 

12 of reproductive cells is always surprisingly high. If somatic cells are very small, 

13 they can outnumber reproductive cells but their biomass is still less than the 

14 biomass of reproductive cells. Only for non-concave benefit functions can the 

15 biomass of somatic cell exceed the biomeiss of reproductive cells. I discuss the 

16 biology of primitive multicellulair organisms with respect to the model predic- 

17 tions. I find good agreement and outline how this work can be used to guide 

18 further quantitative studies of multicelluleirity. 
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1 1 Introduction 



2 Every organism is exposed to mutations that cause variation in inherited traits. Competition 

3 between shghtly different organisms leads to the prohferation of variants that increase fitness. 

4 Most adaptations will fine-tune existing systems but some adaptations lead to new features. 
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undifferentiat ed multicellular organisms can evolve fair l y quickly through mutatio ns of sur- 
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such clonal aggregates do no t have to c ompete against each other for reproduction since 



12 they are genetically identical (iBussl . 



19881 ). This alleviation of reproductive competition has 

13 a profound effect. It allows for a division of labor. Cells can specialize on non-reproductive 

14 (somatic) tasks and peacefully die since their genes are passed on by genetically identical 

15 reproductive cells which benefited from the somatic function. This division of labor turns 

16 multicellular organisms into more than just lumps of cells. They contain cells that are 

17 different in function and appearance. Today a plethora of differentiated organisms exist, 

18 demonstrating the evolutionary success of division of labor. 

19 Most theoretical studies of multicellulari ty analyze the change in level of selection a r id the 



20 consequences 



Michod 



f or reproductive compet i tion 



1997 



Michod fc Roze 



1997 



Busd . 



1988 ; 



Maynard-Smith fc Szathmary 



1997 



200ll ). In this paper I study which conditions make 

22 differentiated multicellularity desirable. When does a differentiated multicellular organism 

23 have higher fitness than an undifferentiated or unicellular organism? In the model presented 

24 here, an organism's fitness depends on a quantitative trait. The quantitative trait is deter- 

25 mined by the number and types of cells in the multicellular organism. The mathematical 

26 model allows one to study which kind of benefits multicellularity must convey to compen- 
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1 sate for its disadvantages. I calculate how much (compared to a reproductive cell) a somatic 

2 cell has to contribute to the quantitative trait to make division of labor advantageous and 

3 determine the optimum number /fraction of somatic cells. 

4 The following section describes the model in detail. In the Results, I will first consider 

5 the evolution of undifferentiated multicellularity. To study the evolution of differentiated 

6 multicellularity I analyze the fitness of organisms of constant size. Thereafter, I study 

7 multicellularity in organisms where the size of the organism and the fraction of somatic cells 

8 is governed by the same evolutionary forces. In the Discussion, I use the insights from my 

9 analysis to discuss a broad spectrum of primitive multicellular organisms. 

10 2 The model 

11 In this work I use the rate of biomass production as a measure of fitness. The rate of biomass 

12 production captures an organism's ability to grow and reproduce. It denotes how much 

13 new biomass per unit of existing biomass an organism can generate per unit of time. For 

14 organisms of equal size, production rates are equivalent to fitness (number of new organisms 

15 produced per organism per unit of time). The model considers how somatic cells, body 

16 size, and benefits of multicellularity affect the rate of production. I distinguish between 

17 reproductive and somatic cells but allow for only one kind of somatic cell. Somatic cells are 
IB different from reproductive cells in that they are terminally differentiated. Their biomass 

19 does not contribute to the next generation. Reproductive cells, on the other hand, contribute 

20 to the next generation. They can be asexually or sexually reproductive. 

21 To derive how somatic cells affect fitness, let us first assume that multicellularity and or- 

22 ganism size have no effect on the rate of production. In this illustrated in Figure [1], a 

23 unicellular organism has the same fitness as a four-cell organism because four unicellular or- 

24 ganisms produce 16 unicellular descendants (16 cells) after two cell divisions, and one four-cell 

25 organism produces four four-cell descendants (16 cells). Indeed, if size and multicellularity 
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1 have no effect on fitness, then all organisms that are entirely composed of reproductive cells 

2 will have the same fitness. Let us use a four-cell organism to derive the cost of somatic cells. 

3 The cost of somatic cells stems from their inability to contribute directly to the next 

4 generation. To quantify this cost, we can compare the rate of production of a four-cell 

5 organism that has two somatic cells with the rate of production of a four-cell organism 

6 without somatic cells. As shown in Figure [H the effect of somatic cells depends on their 

7 size. If the somatic cells are negligibly small, then the rate of production of a differentiated 

8 organism is equal to the rate of production of an undifferentiated organism (forth row in 

9 Fig. [1]). If somatic cells are as large as reproductive cells, then a differentiated four-cell 

10 organism with two somatic cells is able to produce only two new four-cell organisms (8 

11 cells) whereas the undifferentiated organism produces four new organisms (16 cells). Hence 

12 the rate of production of the differentiated organism is 1/2 the rate of production of the 

13 undifferentiated organism. In general, the rate of production (fitness) is reduced by the 

14 fraction of biomass that does not contribute to the next generation. 

In the following we will assign parameters so the biomass and number of somatic and 
reproductive cells to quantify this fitness reduction. I use the Greek letters a and (5 to 
denote the biomass of a reproductive and a somatic cell, respectively, in the adult organism. 
As we will see, it is sufficient to consider the size of a somatic cell relative to the size of a 
reproductive cell. Let B = (3/a denote this size ratio. Let N,. and Ns denote the number of 
reproductive and somatic cells. An adult organism is composed of = A^,. + Ng cells and 
has a body mass of aN^ + (3Ns. Of that biomass (3Ns rests in somatic (sterile) cells and aN^ 
in reproductive cells. Hence, PNs/{aNr + PNg) of the organism's biomass is lost in each 
generation and constitutes the cost of somatic cells. The fitness of an organism with somatic 
cells relative to the fitness of an organism without somatic cells is given by 

1 - (3N,/{aNr + (3Ns) = aNr/{aNr + f3N,) = Nr/{Nr + BN,), (1) 
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ioi B = l3/a as defined above. 

So far I liave assumed tliat tlie rate of production is independent of the organism size 
and that multicellularity does not convey benefits. An overwhelming amount of er npirical 



data shows that the rate of production decreases with the body mass of an organism ([Peters 



19861 ). In particular the annual rate of production per average biomass scales with , 
where W denotes the average body mass of an adult organism and 7 is a scaling factor. This 
relationship holds from unicellular organisms to ma mmals, with b ody masses ranging from 



PetersI (11986 



p. 134), the exponent 7 



approximately 10^^° to 10^ kg. As summarized by 
might range from 0.23 to 0.37. For small organisms 7 is close to 1/4, the typical allometric 
exponent of size. I will therefore use 7 = 1/4 to discuss quantitative results. From above, we 
know that the adult body mass of an organism in this model is given \yy W = aNr + l3Ns. 
Hence, the rate of production decreases by the factor {aNr + f3Ns)~'^. 

I model the advantages of multicellularity as a function of the number of reproductive 
and somatic cells. Let f{Nr,Ns) denote this benefit function. Therefore, the fitness of a 
multicellular organism is given by 

F[Nr, Ng) = cost of somatic cells x cost of size x benefit of multicellularity 

X {aNr + (3Ns)-"' X f{Nr, Ns) 



aNr + PN 
N 

oc 7 '-——f{Nr,Ns), (2) 

where B = (3/a is the size of a somatic cell relative to the size of a reproductive cell. 

This model confirms common sense. If multicellularity does not affect fitness, that is, 
f{Nr,Ns) = constant, then a unicellular organism has a higher fitness than a multicellular 
organism since F{Nr, Ng) < F{Nr, 0) < F{1, 0) for A^^ > 1 and Ng > 0. For undifferentiated 
multicellularity the fitness of an organism is given by N~'^ f{Nr, 0) and multicellularity is 
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only advantageous if it conveys benefits that compensate for the disadvantages caused by 
the size increase. In other words, f{Nr, 0) has to increase more steeply than A^""^ decreases. 

Central to my analysis of multicellularity is the function f{Nr,Ns), which captures the 
benefit of multicellularity. I will assume that somatic and reproductive cells contribute to a 
quantitative trait, x, and that the benefit of multicellularity is a function, f{x), of this trait. 
For simplicity, I illustrate this approach by formulating predator evasion and fiagellation in 
terms of this model. 

For predator evasion, the quantitative trait is given by the size of the organism. It 
increases with the number of cells. Its value determines to which extent the organism is able 
to evade predation. If the organism is big enough, the predator is unable to ingest it and the 
benefit of multicellularity, f{x), is close to 1. For small organisms predation might be severe 
and f{x) close to 0. One can expect a steep increase of f{x) as the organism size surpasses 
the maximum particle size the predator can ingest. Figure [2^ shows a benefit function that 
could be used to describe predator evasion. 

For fiagellation, the quantitative trait is given by the fiagellar drive that the cells of the 
organism provide. The more cells, the more fiagellar drive, which improves the organism's 
ability to maintain its position in a favorable environment. For this example the benefit 
function can be expected to be concave. An initial increase in fiagellar drive might be very 
beneficial by allowing the organism to maintain its position. At some point, however, the 
organism has enough fiagellar drive to maintain its position for most of the time and a 
further increase in fiagellar drive does not yield a substantial benefit. Figure [2b shows a 
benefit function that could be used to model benefits from fiagellation. 

As we will see, my analysis of this model does not require an exact specification of 
the benefit function. But it is necessary to make assumptions about the contribution of 
individual cells to the quantitative trait. 

I am mainly interested in the evolution of early multicellular organisms composed of few 
cells for which it is reasonable to assume (at least as an approximation) that cells contribute 
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additively to the quantitative trait. I use the letters a and b to denote the contribution 
of a reproductive and a somatic cell to the quantitative trait. In an organism with Nr 
reproductive and A^^ somatic cells, the quantitative trait, x, is given by 



X = aNr + bN, 



(3) 



Hence, we have 



f{Nr,N,) = f{aNr + bN,). 



(4) 



1 Let A = a/b denote the contribution of a reproductive cell relative to the contribution of 

2 a somatic cell. Since somatic cells are specialized cells and free of reproductive constraints, 

3 they will generally contribute more to the trait than reproductive cells (6 > a), hence, 

4 A = a/b ^ [0,1]. Furthermore, I consider organisms where somatic cells are usually smaller 

5 than reproductive cells (/? < a) and, hence, -B = e (0, 1]. 

To simplify the analysis, I scale the argument for / so that b = 1, that is, a somatic cell 
contributes one unit to the quantitative trait. By using A = a/b and rescaling the argument 
for /, we can rewrite ([2]) as 



6 Note that if the contribution of reproductive cells to the quantitative trait is insignificant 

7 (i.e., ANr + Ns ^ Ns) then the benefit of multicellularity is simply a function of the number 

8 of somatic cells. 



10 f{x) to monotone increasing and bound functions. Both constraints are reasonable. An 

11 increase of the quantitative trait should not lead to a decrease in fitness, and the fitness of 

12 an organism cannot be increased infinitely. As mentioned earlier, f{x) is scaled such that 

13 lim f{x) = 1. 



F oc 




(5) 



So far I have not made any assumptions about f{x). To simplify the analysis I restrict 



X 



■oo 
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Table [T] summarizes the variables and parameters used in this model. 



2 3 Results 



3 I am interested in the kind of benefit functions, f{x), that promote the evolution of multi- 

4 cellularity. I am also interested in how much a somatic cell has to contribute to the quan- 

5 titative trait (relative to reproductive cells) to compensate for the loss in somatic biomass. 

6 I calculate the optimum fraction of reproductive cells and the optimum fraction of repro- 

7 ductive biomass. This allows us to determine under which conditions undifferentiated or 

8 differentiated multicellularity evolves and to infer on the composition of early differentiated 

9 multicellular organisms. 

10 First, I analyze conditions for the evolution of undifferentiated multicellularity. There- 

11 after I study the evolution of somatic cells in organisms of constant size. I study the uncon- 

12 strained model last. In the unconstrained model the size of the organism and the fraction of 

13 somatic cells is governed by the benefit function f{x). I show that it is possible to calculate 

14 the optimum fraction of reproductive cells and that this fraction is independent of the benefit 

15 function f{x). 

16 3.1 Undifferentiated multicellularity 

The evolution of undifferentiated multicellularity corresponds to an evolutionary transition 
from organisms composed of one (reproductive) cell to organisms composed of several iden- 
tical (also reproductive) cells. Undifferentiated multicellular organisms have, per definition, 
no somatic cells {Ng = 0). For A^^ = 0, the fitness ([2]) simplifies to 




(6) 
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Multicellularity (A^,. > 1) is only advantageous if N''^ f{Nr, 0) > /(I, 0). In other words, the 
benefit function f{Nr,0) has to increase faster than N~'^ decreases. 

I employ this simple case to illustrate the method that will be used to analyze the more 
complex cases. I would like to know for which functions f{x) multicellularity is advantageous 
and what the optimum number of reproductive cells is. To get a general idea of how f{x) 
affects fitness, we can determine the functions f{x) for which the fitness is constant with 
respect to A',.. Let fiso{x) denote these functions. I refer to them as isolines since they 
join points of equal fitness. They are analogous to the lines on topographic maps that join 
points of equal altitude. These isolines can be used to illustrate the fitness landscape with 
respect to benefit functions f{x). The fitness of an undifferentiated organism is given by 
F = N~'^ f (aNr) and constant if f{aNr) oc . Substituting x = aNr, we get 

/iso(a;) oc x'^. (7) 

The gray curves in Figure [3] show these isolines. With the knowledge of these isolines 
it is easy to determine which functions, f{x), promote the evolution of multicellularity. It 
is also simple to determine the optimum number of reproductive cells. An organism with 
Nr reproductive cells has x = aNr as value for the quantitative trait. We have an optimum 
a^opt = ct^r.opt if f{.x) < /iso(a;) for the isoline with /(xopt) = /iso(a^opt)- In a continuous setting 
the optimum satisfies ^/(xopt) = ^/iso(a;opt) for the isoline with f{xopt) = /iso(2;opt)- 

To interpret isoline plots, it might be useful to keep the analogy with topographic maps 
in mind. One can think of x as the distance traveled along a particular trail, f{x), in a 
mountainous region. The highest point, in our case the optimum, Xopt, is reached if the trail 
"brushes" the highest contour line along the trail, /(xopt) = /isola^opt)- None of the points 
along the trail will be above this contour line, f{x) < fiso{x)- 

Figure [3] shows two linear and one concave benefit function. The black bullets indicate 
the optimum for each benefit function. As one can see, a linear benefit function will always 
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1 promote the evolution of multicellularity and will increase f{x) until it reaches a value 

2 close to one. A concave benefit function reaches the optimum earlier and results in smaller 

3 organisms. Undifferentiated multicellularity would not evolve if f{x) increases slower than 

4 the isolines. In particular, multicellularity would not evolve if /(a) > 2^"'. For 7 = 1/4, we 

5 have 2"''' = 0.84 and multicellularity would not evolve if unicellular organisms are able to 

6 benefit from the quantitative trait more than 84% of its full potential. 



7 3.2 Differentiated multicellularity in organisms of constant size 

8 In the previous section I have shown that undifferentiated multicellularity is advantageous 

9 for many benefit functions. In the following I will analyze the evolution of somatic cells 

10 (differentiated multicellularity). For simplicity, I will first analyze the evolution of somatic 

11 cells in organisms of constant size. 

12 From a biological perspective it is relevant to consider organisms of constant size, since 

13 many benefits of undifferentiated multicellularity imply constraints on the size of the organ- 

14 ism. Predator e vasion, for example, is known to promote the evolution of undifferentiated 



15 multicellularity (IBoraas et al. 



19981 ). An organism that uses multicellularity to evade pre- 

16 dation is obviously constrained with respect to size. It has to be larger than the largest 

17 particle that the predator can feed on. Replacing nine large reproductive cells with nine 

18 small somatic cells could decrease its size to dangerous levels. 

19 For simplicity, I will first assume that somatic and reproductive cells have the same size 

20 {B = 1) and explore the more general case of i? < 1 thereafter. 

21 3.2.1 Somatic cells are as large as reproductive cells (B = 1) 

The size of an organism is given by S" = aN,. + (3Ns. If somatic and reproductive cells have 
the same size (5 = /?/« = 1), then the size of the organism can only be held constant if the 
number of cells that compose this organism is constant, that is, N = Nr + Ns = constant. 



11 



In this case the fitness, F{Nr, Ng), depends on one variable instead of two. Using A^^ as this 
variable we can rewrite ([2]) as 



FcxNrf[N-{l-A)Nr], (8) 

because the trait value of an organism with cells is given by x = N — {l — A)Nr. Hence, the 
quantitative trait of an organism of size has a value of at least AN (if A^^ = N) and at most 
A^ - 1 + A (if A^, = 1). The fitness of an organism is constant if / [A^ - (1 - A)Nr] cx A^'^^ 
Expressing the number of reproductive cells, A^^, in terms of the quantitative trait, x, we get 

f\M^{N-x)-\ (9) 



1 Since we rescaled f{x) so that a somatic cell contributes one unit to the quantitative trait, we 

2 can interpret A^ in (l9l) as the value of the quantitative trait of an organism entirely composed 

3 of somatic cells. Hence, N — x could be interpreted as the value by which the quantitative 

4 trait is reduced due to the existence of reproductive cells, i.e., the cost of reproductive cells 

5 in terms of the quantitative trait. 

Let us first calculate the optimum number of reproductive cells for a linear benefit func- 
tion, f{x) = cx. To guarantee that < /(x) < 1 for all organisms composed of A^ cells, we 
constrain c to < c < ^ 1/N. As mentioned above, the optimum, Xopt, satisfies 



^f{xopt) = ^/iso(xopt) for the isoline with /iso(a;opt) = f{xopt)- Hence, we have to solve 



the two equations ^/(xopt) = ^/iso(a;opt) and /iso(a;opt) = /(a^opt) for Xopt and the (irrele- 
vant) constant k in fiso{x) = k{N — x)^^. Solving these equations for the benefit function 
f{x) = cx, we get Xopt = N/2 which corresponds to and optimum number of reproductive 
cells of 

As we can see, the optimum value for the quantitative trait is independent of A. The 
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1 parameter A does, however, determine how many somatic cells are necessary to reach the 

2 optimum value Xopt for the quantitative trait, and hence A^r-.opt- Remarkably, the optimum 

3 number of reproductive cells is independent of the slope, c, of the linear benefit function. 

4 We only constrained the slope so that the linear function does not exceed one {f{x) < 1). 

5 Equation ffTOj) also shows that the number of reproductive cells is usually greater than N/2 

6 and only equal to N/2 if reproductive cells do not contribute to the quantitative trait {A = 0). 

7 If the optimum number of reproductive cells, Nj-^opt, is less than N, then somatic cells are 

8 advantageous and differentiated multicellularity is likely to evolve. From ( |T0|) we see that 

9 this is only the case if A < 1/2 which means that somatic cells have to contribute twice as 

10 much as reproductive cells to the quantitative trait to justify their existence. 

11 We can summarize our results for organisms of constant size, uniform cell sizes, and 

12 linear benefit functions: (a) such organisms contain many reproductive cells, and (b) somatic 

13 cells in such organisms have to contribute substantially more to the quantitative trait than 

14 reproductive cells. In the following I analyze isoline plots to illustrate that this result holds 

15 for many nonlinear benefit functions. 

16 Figure UK shows four benefit functions and isolines fiso{x) = [N — x)~^ for = 32. 

17 The solid line represents a linear benefit function that satisfies the requirements from above 
IB (/(x) < 1). It is evident that the isolines and the benefit function have the same slope 

19 at Xopt = N/2 = 16. The Figure illustrates that Xopt does not depend on A (the relative 

20 contribution of a reproductive cell to the quantitative trait). Even though Xopt is always 

21 with respect to A, this parameter determines how many somatic cells, if any, are required 

22 to reach this Xopt- The top axis shows how x corresponds to A^^ for four different values of 

23 A. As A increases to one, the range of possible values for x, AN to N — {1 — A), shrinks (to 

24 the right). This is no surprise. If the contribution of reproductive and somatic cells to the 

25 quantitative trait are about the same, then the total value of the quantitative trait changes 

26 little if one substitutes a reproductive cell for a somatic cell. One can see that if A is larger 

27 than 1/2, then the quantitative trait of an undifferentiated organism already exceeds the 
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optimum value of Xopt = N/2 < AN for a linear benefit function. There is no need for the 
organism to evolve somatic cells. 

Figure nil contains two concave benefit functions (dashed and dash-dotted curves). It is 
easy to localize Xopt for these functions and obvious that their Xopt is smaller than the Xopt 
for the linear benefit function. For a concave benefit function, organisms need less somatic 
cells to optimize fitness and the functional demand on somatic cells to justify their existence 
increases {A has to be even smaller, see upper axis in Fig. H^i). Hence, the results from above, 
that organisms have few somatic cells and that somatic cells have to contribute substantially 
more to the quantitative trait does also hold for concave functions. 

Examining the isolines in Figure H^i we see that only a convex function can lead to 
organisms with many somatic cells and few reproductive cells. A convex function would 
describe a situation in which the organism has to obtain a minimum threshold value to 
benefit from the quantitative trait. In such a situation the functional demand on somatic 
cells is relaxed and organisms might require many somatic cells to optimize fitness. 

3.2.2 Somatic cells that are smaller than reproductive cells {B < 1) 

If somatic cells are smaller than reproductive cells, then the total number of cells, A^, can 
change even if the size of the organism remains constant. If, for example, reproductive cells 
are twice as large as somatic cells {B = 1/2), then one reproductive cell can be replaced by 
two somatic cells, which increases the total number of cells by one but keeps the organism 
size constant. 

If somatic cells are half the size of reproductive cells, we can simply apply the results 
from above by changing the parameter A. Since contributions are additive in my model, 
somatic cells of size (3 that contribute b to the quantitative trait are equivalent to somatic 
cells of size f3/2 that contribute 6/2 to the quantitative trait. In other words, if somatic cells 
that are as large as reproductive cells are beneficial, then somatic cells that are half as large 
and contribute half as much to the trait have to be beneficial as well. Hence, to get result 
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for S < 1, we only have to consider the results from above and replace A with AB. For 
example, for a linear benefit function the optimum number of reproductive cells is given by 



1 and somatic cells are advantageous if AB < 1/2. For concave benefit functions, we can 

2 conclude that more than half of the biomass of the organism will rest in reproductive cells 

3 and that somatic cells need to satisfy AB < 1/2 to justify their existence. Appendix A 

4 contains a more technical and detailed analysis of the case B < 1. For completeness, I 

5 analyze a model in which the number of cells is held constant (as opposed to the size of the 

6 organism) in Appendix B. 

7 3.3 The complete (unconstrained) model 

8 Let us now study the unconstrained model. By not restricting the number of cells or the 

9 size of the organism, I assume that its size and the optimum fraction of reproductive cells 

10 are governed by one evolutionary force. In this case the quantitative trait, x, can no longer 

11 be expressed as a function of A^^- It depends on Nj. and A^^. This makes the calculation and 

12 visualization of isolines unwieldy. Instead, we can actually calculate the maximum of the 

13 fitness function. In the following we will realize that If the size and the composition of the 

14 organism can change freely, then the optimum fraction of reproductive cells is independent 

15 of f{x) and can be calculated. 

Using Ng — N — Nj. and q — Nj./N we have 




(12) 
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and can rewrite ([2]) as 



^ ^ ^nr. w^ T^^f[N{l-{^- A)q)\. (13) 

In the following I assume that there is at least one reproductive cell (g > 0) and that somatic 
cells contribute more to the quantitative trait than reproductive cells {A < 1). Further, I 
assume that /(x) is differentiable and monotone increasing > 0). We can calculate 

df df dx df 

and 

dq dx dq dx 

Since none of the factors in f |T^ and ( |T5l) equal zero, we can express in terms of 

df df Nil -A) 



dq dNl-{l-A)q 
Applying the product rule of differentiation to (fT3l) we get 



(16) 



dF if df 

m'^-J^^dN- ^'^^ 

If the number of cells, A^, is optimal, then |^ = and hence §^ = which can be 
substituted into ( |T6i) to give 

^ " ~l-{l-A)q^- ^^^^ 
Differentiating ( fT3l) with respect to q results in 

dF d f 

— (X [(1 - B)q + B]f - (1 + 7)g(l - B)f + g[(l - B)q + B]^ (19) 
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Using (fTS!) we can substitute ^ in (fT9|) and get 
dF 



dq 



/■ (20) 



The fraction of reproductive cells is optimal if ^ = 0. Since / > 0, the optimum 
fraction of reproductive cells can be determined by calculating for which q the factor in (l20l) 
equals zero. Multiplying this equation by 1 — (1 — A)q shows that terms quadratic in q, i.e., 
{1 — A){1 — B)q^ cancel. Hence, this equation is linear in q and can be solved to give 

Thus, we are able to calculate the optimum fraction of reproductive cells, gopt- Remark- 
ably, gopt is independent of the benefit function /(x). The benefit function will, however, 
determine the size of the organism. 

Figure [5] shows an implicit plot of gopt as a function of the parameters A and B. The 
curve for gopt = 1 gives the threshold for parameter values that favor the evolution of 
differentiated multicellularity. This curve is given by AB = (= 1/5 for 7 = 1/4). In 
particular, if somatic cells are as large as reproductive cells {B = 1), then they are only 
beneficial if A < For 7 = 1/4, somatic cells have to contribute five times as much to the 
quantitative trait than the reproductive cell. Similarly, if somatic cells contribute as much 
to the quantitative trait as reproductive cells {A = 1), then, to be advantageous, their size 
has to be a fraction of the size of reproductive cells. In particular, this fraction has to be 
less than (1 + 7)/7. 

Figure [5] shows that for small somatic cells (small B) , the ability of reproductive cells to 
contribute to the quantitative trait {A) has little effect on gopt- For small B, the denominator 
in (I2TI) is approximately 1 and gopt ~ 7^^-B. From (fT3!) we also see that A appears in the 
equation for F only in the term 1 — (1 — A)q. If g is small (because of small B), then 
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1 1 — {1 — A)q ^ 1 and A has little effect on the fitness of the organism. 

2 That the effect of A (the ability of reproductive cells to contribute to the quantitative 

3 trait) on F and ^r.opt depends on B (the size of somatic cells) is an important result for 

4 our understanding of the evolution of somatic cells. Many reproductive cells have to grow 

5 to a minimum size before they can initiate cell division. Newly evolved somatic cells are 

6 presumably as large as reproductive cells but are instantaneously relieved of reproductive 

7 size constraints. The minimum size at which a somatic cell can still function might be much 

8 smaller than the minimum size of a reproductive cell. Organisms will have the tendency to 

9 evolve somatic cells that are as small as possible (decrease B). Equation [2T] shows that a 

10 decrease of B increases the optimum number of somatic cells (decreases qr,opt)- This will 

11 further increase the selective pressure to reduce the size of somatic cells because the organism 

12 has now more somatic cells that should not be unnecessarily large. This feedback loop might 

13 continue until there are many, small somatic cells. At this point (small B) the contribution 

14 of reproductive cells to the quantitative trait has no major effect on fitness (see Fig. |8b)- 

15 Hence, reproductive cells can cease to contribute to the somatic function (the quantitative 

16 trait in my model) with little effect on fitness. They are free to dedicate their existence 

17 fully to reproductive duties. Such an evolutionary feedback loop promotes the evolution of 
IB organisms with a strict division of labor between many, small somatic cells and few, large 

19 reproductive cells. 

20 It is important to emphasize that we treated q and as continuous variables. Especially 

21 for small multicellular organisms they are, however, discrete. For example, a bicellular 

22 organism can have a g of 1/2 or 1. From Figure [5] and (pTj) we know that somatic cells in 

23 such a bicellular organism are (even for A = 0) only advantageous if they are much smaller 

24 than reproductive cells (5 < 7). Also, for B = 1 the optimum fraction of reproductive cells 

25 is always larger than 1/(1 + 7) (= 4/5 for 7 = 1/4). Hence, evolutionary transitions to 

26 differentiated multicellularity with somatic cells that are as large as reproductive cells can 

27 only happen in organisms that are composed of at least six cells. 
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4 Discussion 



I have presented a model for the evolution of undifferentiated and differentiated multicellu- 
lar organisms. In my model three factors determine an organism's fitness: (a) its size (or 
biomass), (b) its investment in somatic (terminally differentiated) cells, and (c) a quanti- 
tative trait that is mainly determined by the number and kind of cells that the organism 
is composed of. The quantitative trait, x, affects the fitness of the organism via a benefit 
function, f{x) (see Fig. [T]). For simplicity I assume that the cells of a multicellular organism 
contribute additively to the quantitative trait. Since somatic cells are specialized and ter- 
minally differentiated, they can contribute more to the quantitative trait than reproductive 
cells. 

I analyze under which conditions (benefit function, contribution to the quantitative trait, 
size of somatic cells, etc.) the evolution of undifferentiated and differentiated organisms 
is favored, and calculate the optimum fraction of somatic cells. My analysis shows that 
undifferentiated multicellularity is favored by many benefit functions. The evolution of 
undifferentiated multicellularity is, however, unlikely if the unicellular organism is already 
able to receive large benefits from the quantitative trait. In particular, multicellularity will 
not evolve if the unicellular organism benefits from the quantitative trait more than 84% 
(= 2-^ for 7= 1/4) of its full potential. 

My model suggests that primitive differentiated organisms will generally have a small 
fraction of somatic biomass. If somatic cells are as large as reproductive cells and the benefit 
function linear or concave, then the fraction of somatic cells is always less than or equal to 
1/2. If somatic cells are smaller than reproductive cells, they might occur in large numbers, 
but their biomass will still be at most 1/2 of the total biomass. Somatic cells compose more 
than 1/2 of the organism only if the benefit function is convex. 

In the following I discuss the biology of primitive multicellular organisms. First I discuss 
undifferentiated, then differentiated organism. I use experimental data from volvocine algae 
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to demonstrate how experimental observations can be compared with model predictions from 
this work. At the end I point out the limitations of my model. 

Most benefits of undifferentiated multicellularity relate to an organism's ability to evade 
predators or its ability to secure a favorable position in the environment. Predator evasion is 



5 comrnonly recognized as a driving force 



(IBuss , 



1988 



King 



20M). 



Boraas et al. 



br th e evolution of undifferentiated multicellularity 



(119981 ) showed that unicellular algae can evolve mul- 



ticellularity within few generations after exposure to a phagotrophic predator. Phagotrophic 
and many other predators face an upper size limit for the particles they can ingest. A simple 
"sticking together" of cells provides protection by exceeding these size limits. In this case 
the quantitative trait is the size of the organism. If it exceeds a certain value, the organism 
benefits substantially from it (see Fig. [2^). 

Multicellularity is also known to improve an organism's ability to obtain a favorable po- 
sition in the environment. In particular, the flagellation constrain t dilem ma is believed to 



play an important role in the evolution of multicellularity ( IMargulis 



19811 1 . Many eukaryotic 



15 cells face the dilemm a that they are unable to r naint a in flagellation during cell division and 



16 hence, lose motility (iBonnerl . 



1965 



Margulis 



1981 



Buss 



1988 



Koufopanou 



1994 



Kirk 



19971 ). In an undifferentiated multicellular organism, motility can be maintained. Multicel- 



18 lularity can a.. 



so increase the speed of a n organism. Many cells can provide more drive than 



19 a single one (jSommer fc Gliwicd . 



19861 ). For this example, the quantitative trait, flagellar 
drive, determines the organism's ability to reach a favorable position in the environment 
which constitutes a benefit (see Fig. [2b). 

Multicellularity can also improve an organism's ability to float. Many algae lack flag- 
ella. They regulate buoyancy throu gh the production of carbohydrate ballast and/or gas 



24 inclusions ( 1 Graham fc Wilcox 



19991 ). Filamentous growth in combination with the secre- 
tion of extracellular polymeric substances allows the formation of mats that provide a sta- 



26 ble structure which can be used to regulate buoyancy by trapping bubbles (jPhillipd . 



Graham &: Wilcox 



1958 



19991 ). In this case the quantitative trait might be given by the tightness 
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of the mat. Tight mats allow to trap many bubbles and allow the cells in that mat to stay 
close to the surface water where they receive more light. 

Let us now consider differentiated multicellularity. My analysis predicts that primitive 
differentiated multicellular organisms will generally have many reproductive cells. More 
precisely, in most cases I would expect more than 1/2 of an organism's biomass to rest in 
reproductive cells. In the following I discuss algae and slime molds, two groups of organisms 
for which quantitative data exist. 

Volvocine algae are an excellent group of organisms to study differentiated multicellu- 
larity. The ir mu l ticellu lar complexity ranges from undifferentiated to highly differentiated 



organisms (IKirk 



19971 ). The most primitive differentiated forms have somatic cells that 

11 maintain flagellation during cell division of reproductive cells. The flagellar beating is also 

12 important to provide a constant nutrient supply. It stirs the medium and prevents a nu- 

13 trient depletion of the organi sm's boundary lay ers which would occur due to the nutrient 



upta ke by the o rganism itself (jSolari et al. 



15 esi s (iBell 



by 



Solari et al. 



20061 ). According to the source-and-sink hypoth- 



1985 ) somatic cells can also increase the uptake rate of nutrients, but experiments 



(l2006l ) suggest that the stirring of the medium plays a more important role 
in nutrient supply. 

The smallest di f ferentiated co l onies in volvocine algae have 32 cells and 4, 8, or 16 are 



19 somatic ( IGoldstein 



1967 



Bonner 



2003bl ). This is in nice agreement with the model predic- 



tion. Allometric data about soma and ger m in volvocine alga e shows that there is more germ 



21 tissue than somatic tissu e in all species f|Koufopanou 



22 use the data collected by 



1994J . Fig. 7). Furthermore, we can 



Koufopanoul (Il994j ) to calculate q and B. In Figure [6l I compare 



this data with the model predictions (12T]) . The good agreement between the model and the 
data suggests that the size and fraction of somatic cells in volvocine algae are governed by 
the same benefit function. 

To compare the model predictions with the experimental data, I assumed that reproduc- 
tive cells do not contribute to the quantitative trait {A = 0). In this case the benefit of 
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multicellularity is modeled as a (largely arbitrary) function of the number of somatic cells. 
To calculate the optimum fraction of somatic cells, it was not necessary to specify this func- 
tion. My model is in this case surprisingly general and should be applicable to a wide range 
of primitive differentiated multicellular organisms. 

Figure [6] suggests that the allometric exponent for volvocine algae, 1/5, is smaller than 
the typical allometric exponent of 1/4. It should be possible to determine experimentally 
the allometric exponent of volvocine algae and compare my predi ction with e x perim ental 



Solari et al. 



mm . one 



observations. From experiments similar to the ones conducted by 
could also learn something about the shape of the benefit function by manipulating the 
number of (functional) somatic cells. 

Another organism group that is commonly used to study primitive differentiated multicel- 
lularity is the slime molds. Slime molds such as Dictyostelium discoideum feed as individual 
cells until food becomes scarce, at which point they form a multicellular mass that migrates 
to a suitable spot and differentiates into a fruiting body. The somatic stalk of the fruiting 



15 body lifts the spores above the ground to facilitate more efficient dispersal (IBonnerl . 



19671 . 



2003al ). In this case, the quantitative trait that conveys the benefit of multicellularity is 



the height of the stalk. The higher the stalk the more efficient is the dispersal. According 
to my model, we would expe ct the biomass of t he stalk to be less then 50% of the total 



19 biomass of the fruiting body. 



FarnsworthI (119751 ) measured the percentage (dry weight) of 



stalk as a function of the temperature during culmination. This percentage changes from 
about 20% at 18°C to 13% at 27°C. Hence, most of the fruiting body is indeed composed 
of reproductive cells. The data suggest that somatic cells are slightly more advantageous 
at lower temperatures since the fraction of somatic cells increases. Similarly, in Myxococcus 



24 xanthus, a 



ruiting body forming 



are spores (lO' Connor fc Zusman 



jacte rium, more than 61% of the cells in a fruiting body 



1991 



Table 2). 

It would be interesting to collect data about soma and germ in slime molds that is 
analogous to the data collected for volvocine algae. A comparison of such data with the 
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1 presented model would be of particular interest since the quantitative trait is most likely 

2 the size of the stalk. This trait is easy to measure and would allow conclusions about the 

3 benefit function. The model presented in this p aper should g uide the researcher in their 



Koufopanod (Il994l ) reported the average 



4 data collection and presentation. For example, 

5 and standard deviation of the number and size of somatic and reproductive cells. In the 

6 light of my analysis it seems to be of greater biological importance to report the average 

7 and standard deviation of the fraction of reproductive cells, g, and the size of somatic cells 

8 relative to reproductive cells, B. 

9 For most conspicuous organism such as plants and animals the number and biomass of 

10 somatic cells vastly outnumbers that of reproductive cells. Notably, all of these organisms 

11 are much more complex than the primitive multicellular organisms that are the focus of 

12 this study. They contain many somatic cell types that form organs and interact with each 

13 other in complex ways. It is important to keep in mind the kind of organisms that the 

14 model is able to describe. I make two key assumptions: (a) The benefit of multicellularity 

15 can be modeled as a function of a quantitative trait. In particular, for a given value of the 

16 quantitative trait, the benefit does not depend on the number of reproductive cells, and (b) 

17 that cells contribute additively to that trait. If one of these two assumptions is not satisfied, 
IB the results can be quite different. 

19 In Appendix C I analyze a model in which the fitness of an organism depends on how 

20 much of a limiting resource (e.g.. Nitrate) each reproductive cell obtains. This resource, 

21 after it has been acquired by the (maybe multicellular) organism, has to be divided between 

22 reproductive cells. This is in disagreement with assumption (a) since the benefit depends 

23 on the number of reproductive cells (the less reproductive cells, the more nutrients each 

24 reproductive cell receives). As shown in Appendix C, such a situation does not favor the 

25 evolution of undifferentiated multicellularity. Differentiated organisms will tend to be small 

26 and are mostly composed of somatic cells. 

27 One might also wonder how my model can be applied to organisms for which the distinc- 
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1 tion between somatic and reproductive cells is not so clear cut. For primitive differentiated 

2 multicellular organisms, "somatic" cells can be characterized by a delayed cell division or a 

3 reduced probability of reproduction, rather than no cell division or no reproduction at all. It 

4 is straight forward to incorporate this developmental plasticity into my model by modifying 

5 the term that captures the cost of somatic cells. This cost is given by the biomass that is 

6 lost due to the existence of somatic cells (or more generally the resources that are lost) . For 

7 terminally differentiated cells this is just given by the biomass of the somatic cells but can 

8 be modified to reflect any developmental plasticity. If, for example, "somatic" cells have 

9 approximately a 50% chance of reproduction, then the average evolutionary cost of somatic 

10 cells is given by 50% of the somatic biomass. 

11 In this work I mathematically described the costs and benefits of differentiated and 

12 undifferentiated multicellularity. I showed that multicellularity can evolve readily if cells of a 

13 multicellular organism contribute additively to a quantitative trait that benefits the organism 

14 in a manner that is independent of the number of reproductive cells. Multicellularity is 

15 especially beneficial if a single-cell organism alone cannot benefit from the quantitative trait 

16 substantially. Only if the single-cell organism is able to exploit the quantitative trait to 84% 

17 (= for 7= 1/4) of the quantitative traits full potential will multicellularity not evolve. 

18 I showed that evolutionary forces that are based on such quantitative traits will generally 

19 evolve multicellular organism with few somatic cells even if somatic cells contribute much 

20 more to the quantitative trait than reproductive cells. 

21 In particular, for the complete model (organism size and fraction of somatic cells is 

22 determined by the benefit of multicellularity) and for somatic cells that are as large as 

23 reproductive cells, the optimum fraction of somatic cells is always less than 7/(1 + 7) = 

24 1/5. As a consequence, under such conditions multicellular organisms can only benefit from 

25 somatic cells if they are composed of at least five cells. Somatic cells can be numerous if 

26 they are very small compared to reproductive cells but their biomass will still be less than 

27 that of the reproductive cells. In the presence of many, small somatic cells, the contribution 

24 



of reproductive cells to the quantitative trait has httle effect on the fitness of the organism. 
This allows reproductive cells to specialize on the reproductive function and paves the way 
for a strict division of labor between reproductive and somatic cells. 
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1 Appendix A: Optimum number of reproductive cells in 

2 organisms of constant size with small somatic cells. 

3 In this section I derive the results for organisms of constant size S and somatic cell that are 

4 smaller than reproductive cells {B < 1). The unit for S is chosen so that an undifferentiated 

5 organism is composed of S reproductive cells. The quantitative trait of such an undifferen- 
5 tiated organism totals AS. A differentiated organism of constant size with one reproductive 

7 cell has {S — 1)/B somatic cells and its quantitative trait equals {S — 1)/B + A. Notably, an 

8 organism's quantitative trait can range from AS to {S — 1)/B + A and depends on A and 

9 B. 

For constant size the fitness ([5]) is given by 

F oc Nrf{ANr + Ns) (22) 
= Nrf[S/B-{l/B-A)Nr], (23) 



where we used A''^ = {S — Nr)/B. Isohnes are given by 

fM(xiS/B-x)-' (24) 



10 and are independent of A. The term S/B can be interpreted as value of the quantitative 

11 trait that an organism entirely composed of somatic cells would have. 

For a linear benefit function f{x) = cx we can calculate the optimum as Xopt = S/{2B) 
which corresponds to 

12 Reproductive cells will constitute 1/[2(1 — AB)] > 1/2 of the biomass of the organism. As 

13 for B = 1, most of the organism's biomass will be reproductive cells. Somatic cells are only 

14 beneficial if AB < 1/2. 
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Figure [7] shows the isohne landscape and three benefit functions. Figure [7h illustrates 
how the isolines depend on the parameter A and B. Figure [7^ shows how the correspondence 
between x and A^^ depends on A and B. For example, an increase of A would move the upper 
left end of the line that maps x on A^^ to the right and steepen its slope. Increasing B would 
move the lower right point of this line to the left and also steepen the slope. It would also 
change the isolines which approach S/B asymptotically. As for constant A^ most benefit 
functions and parameter combinations will lead to a fairly large number of reproductive cells 
or a large fraction of reproductive biomass. Only for convex benefit functions would A^r.opt 
be small. 



Appendix B: Optimum number of reproductive cells in 
organisms with a constant number of cells. 

In this section I analyze the evolution of small somatic cells (i? < 1) in organisms that are 
composed of a constant number of cells (A^ = constant). The fitness ([5]) is given by 

F oc -T—f fA^ - (1 - A)Nr] , (26) 

[Nr + B{N - Nr)Y+^' ^ ^ ^ '^J' ^ ^ 

and constant if / [A^ - (1 - A)Nr] oc N-^[l + B{N/Nr - 1)]^+^. Expressing A^^ in terms of 
X, we get 

IS — X 

We can interpret A^ — x as the decrease of the quantitative trait due to the existence of 
reproductive cells, and x — AN as the increase in the quantitative trait (compared to undif- 
ferentiated organisms) due to somatic cells. 

Let us now study how a change in the size of somatic cells affects the optimum number 
of reproductive cells. The isolines are given by 1/{N — x) for B = 1. For i? < 1, we can 
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rewrite equation (1271) as 



[iV 



1-BA 
1-B 



1+7 



/iso(a;) oc 



N -X 



(28) 



1 and notice that fiso{x) approaches (A^ — x)'^ for S — > 0. The shape of the isohne is entirely 

2 determined by the factor {1 — AB)/{1 — B) and different combinations of parameters A and i? 

3 can result in the same isoline. Since A appears only in the term 1 — AB, it has less influence 

4 on the shape of fiso{x) if B is small. This can be explained intuitively. If somatic cells are 

5 very small, they are not very costly and how efficient they are (compared to reproductive 

6 cells) is less important. 

7 Figure (Hti shows isolines for different parameter combinations. I choose A = 0.1, 0.25, 

8 0.5, and 0.75, and B = 0.02, 0.1, 0.2, 0.5, and 1. It illustrates the analytical results. 

9 Isolines change from (A^ — x)^^ (purple line) to (A^ — x)'' (gray line) and different parameter 

10 combinations can result in similar isolines. The parameter A has little affect on the isoline 

11 if somatic cells are small (small B). The top axes of Figure [8^ show how x corresponds 

12 to Nr for different values of A. Interestingly, if B is small enough, then isolines can have a 

13 negative slope. In other words, even a constant benefit function would promote the evolution 

14 of somatic cells. We can calculate that the slope of fiso{x) is negative at xpf^=N = 

15 if i? < 7/(1 + 7) = 1/5. For somatic cells of that size the disadvantage of loosing a 

16 reproductive cell is compensated for by the size decrease (smaller organisms have higher 

17 rates of production). We will encounter this threshold again during our analysis of the 

18 complete (unconstrained) model. 

19 Figure [HJo shows the fitness, F, as a function of N^. I plot F for A = 0.1 and indicate 

20 the maxima for the other A values. As expected, a decrease in B leads to an increase of F 

21 and a decrease of A^r-,opt- Also, F and A^r-.opt increase with A, but less so if B is small. 
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Appendix C: Multicellularity in organisms in which the 
benefit of multicellularity depends on the number of re- 
productive cells. 

In this section I analyze a small but significant variation of my model. To model the benefit 
of multicellularity, I assumed that this benefit is a function of a quantitative trait to which so- 
matic and reproductive cells contribute additively. In this section I analyze a model in which 
the benefit of multicellularity (for a given value of the quantitative trait) does also depend on 
the number of reproductive cells. In this version reproductive and somatic cells contribute 
(additively) to the acquisition of a resource that is desperately needed by reproductive cells. 
The more a reproductive cells has of this resource, the faster it can grow and the larger is 
the probability of survival and, consequently, the fitness of the organism. Let f{x) denote 
the benefit from this resource if the organism manages to supply each reproductive cell with 
an amount x of the resource. 

If one somatic (reproductive) cell acquires j3 (a) of the resource, then a total of aNr+jSNg 
can be allocated between the reproductive cells and each cell would receive The 
benefit of multicellularity is then given by 



In the following, I analyze this fitness function. I show that such benefits of multicel- 
lularity are (a) an unlikely source for the evolution of undifferentiated multicellularity, (b) 
differentiated organisms of constant size would have few reproductive cells, and (c) if organ- 




(29) 



and the fitness of the organism given by 




(30) 
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1 ism size and the fraction of somatic cells are governed by f{x) then the optimum number of 

2 reproductive cells is given by one. 



3 Undifferentiated Multicellularity 

4 It is easy to see that for A^^ = 0, equation fl30l) simplifies to F = N~'^f{a). Fitness, F, 

5 is strictly monotone decreasing with A^^ and optimal for A^^ = 1- Obviously, if all cells 

6 contributed linearly to the acquisition of a resource and divide that resource equally among 

7 each other, then each cell gets as much as it would get if it were on its own. In other words, if 

8 cells contribute linearly to the acquisition of a resource that is apportioned among them, then 

9 multicellularity conveys no advantages and a unicellular organism has the highest fitness. 

10 Limited resources that have to be divided between reproductive cells can only trigger the 

11 evolution of multicellularity if cells act synergistically. The multicellular organism has to be 

12 more than just the sum of its parts. 

13 Constant size 

For constant size S = aN^ + PNg, the fitness fl30|) is given by 



The amount of resources that each reproductive cell receives depends on A and B. It can 
range from X]y^=s = ^ to X]y^=i = A + {S — 1)/B. The isolines are given by 



F cx Nrf{A + Ns/Nr) 



(31) 



Nrf [S/{BNr) - (1/5 - A)] . 



(32) 



/iso(a;) (xx + l/B-A. 



(33) 



14 As one can see, the isolines are linear. For linear benefit functions f{x) = cx with f{S/ B) < 1 

15 the optimum number of reproductive cells equals one if AB < 1. For B = 1, this simplifies 
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to A < 1. More so, once differentiated multicellularity is advantageous the optimum number 
of reproductive cells is given by one. Figure [9] shows isolines and the non-linear relation 
between x and Nj. for A = 0.1. I plot three possible benefit functions. It is obvious that only 
very steeply increasing benefit functions would lead to organisms with many reproductive 
cells. Most benefit functions would result in few reproductive cells. 

The complete Model 

Let us now analyze the unconstrained model. As in the main text, the optimum size of the 
organism as well as the optimum number of reproductive cells are governed by the benefit 
function f{x). As we will see, we can use this to our advantage and calculate the optimum 
number of reproductive cells analytically. We can even calculate the isolines with respect to 
size and conclude that benefits from resources that have to be allocated between reproductive 
cells are an unlikely cause of the evolution of multicellularity. 
The fitness is given by 




(34) 



If there is an optimum, it has to satisfy = 0. We have 



dF 



Nr -(1 + ^)B{Nr + 5Ar5)-(2+7)/ + (AT^ + BNs)-''^+"'^N; 



— 1 




(35) 



and hence 




(l+7)i?(iVr + 5Ars)-(2+7)/ 



(36) 
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at optimum, which can be used to simphfy 
OF 



dNr 



-N,{Nr + BNs)-^'+^^N-'^ (37) 



to 



^ = - [(1 + 7)(A^. + BNs)' - 1] {Nr + BNsr^'+'^^f (38) 

which is always negative. Hence, whenever the fitness is constant with respect to changes in 
Ng, Nr will decrease. Since there can never be less than one reproductive cell the optimum 
number of reproductive cells is given by iV^^opt = 1- The fitness of the corresponding organism 
is given by 

F oc (1 + B{N - 1))-M/[A^ - (1 - A)]. (39) 
The isohnes with respect to N are given by 

Uo{^) = [x + l/B-A]'^\ (40) 



1 Hence, multicellularity based on divisible resources is only advantageous for convex benefit 

2 functions that grow faster than x^'^'^. The resulting organism would contain one reproduc- 

3 tive cell and would generally be fairly small (contain few cells). This suggests that benefits 

4 from resources that have to be allocated between reproductive cells are an unlikely cause of 

5 the evolution of multicellularity. 
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Figure 1: Cost of somatic cells in differentiated organisms. If the rate of production is 
independent of organism size, then every undifferentiated organism has the same fitness. 
Four unicellular organisms produce 16 cells after two cell division (first row), as does one 
four-cell organism (second row). A four-cell organism with two somatic cells (third row) 
produces half the biomass of an undifferentiated organism, because only half its cells are 
able to contribute to the next generation. If the two somatic cells are extremely small 
(fourth row) then nearly all the biomass is concentrated in reproductive cells and fitness is 
equal to the fitness of an undifferentiated organism. In general the fitness of a differentiated 
organism relative to an undifferentiated is given by the fraction of biomass that is used for 
reproduction, that is, by Nr/{Nr + BNg) = 1 — BNs/{Nr + BNs), where B is the size of a 
somatic cell relative to the size of a reproductive cell. 
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Figure 2: The use of benefit functions to model the benefit of multicellularity. (a) Predator 
evasion: Multicellularity can be beneficial by allowing the multicellular organism to evade 
predation. We can model this benefit as a function of the organism size (a quantitative 
trait). Multicellularity is not beneficial if the size of the organism does not exceed the 
predators upper prey size limit. Large benefits can be expected for size increases that 
surpasses this prey size limit. Once this threshold is exceeded and the organism immune 
against predation, further size increases bring httle additional benefits, (b) Flagellation: 
Multicellularity can also increase the motility of an organism by increasing its fiagellar drive 
(a quantitative trait). Initial increases in motility can bring large benefits because they 
increase an organism's ability to reach favorable environments. At some point, however, a 
further increase in fiagellar drive brings little benefits because the organism is already able 
to access the most favorable environment. 
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Figure 3: Optimum number of (reproductive) cells in undifferentiated organisms. The gray 
curves show isolines, i.e., benefit functions, fiso{x) = x~'^ (for 7 = 1/4), for which the 
fitness ([5]) of the organism is constant with respect to changes in the number of cells (see 
Eq. [7]). The solid and dotted lines show linear benefit functions. The dashed curve shows a 
concave benefit function. The number of reproductive cells is optimal if f{x) < fiso{x) for 
the isoline with /(xopt) = /iso(2^opt)5 where Xopt = a-^r,opt- Black bullets indicate the optima, 
(^^opt, /iso(a^opt)), for the given benefit functions. 
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Figure 4: Optimum number of reproductive cells in organisms of constant size (A^ = 32) 
with somatic cells that are as large as reproductive cells {B = 1). (a) The gray curves show 
isolines with respect to changes in A^^ (see Eq. |9]). Isolines are independent of A for B = 1. 
But, since x = N — {1 — A)Nr, the parameter A affects the relation between A''^ and x. 
The top axis illustrates how this relation changes with A. (b) Fitness of an organism as a 
function of A^^- The curves correspond to the benefit functions from above for A = 0.1. To 
keep the figure uncluttered, I do not plot the fitness curves for A = 0.25, 0.5, and 0.75, but 
indicate only the maxima. 
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Figure 5: Optimum fraction of reproductive cells. If the size of an organism and the fraction 
of reproductive cells, q, can adjust freely, then the optimum fraction of reproductive cells is 
independent of f{x) and given by (I2T]) . The optimum fraction depends on A (the contribution 
of a reproductive cell to the quantitative trait relative to the contribution of a somatic 
cell), B (the size of a somatic cell relative to the size of a reproductive cell), and 7 (the 
allometric exponent). This implicit plot shows (for 7 = 1/4) which values A and B result in 
Qopt = 1, 0.75, 0.5, 0.25, and 0.125. Differentiated multicellularity for an organism of size N 
is beneficial if gopt < 1 — 1/A^. Note that A becomes irrelevant for small B. 
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Figure 6: Fraction and size of somatic cells in volvocine algae. I plot the fraction of repro- 
ductive cells, q = Nr/N, as a function of the size of somatic cells relative to reproductive 



Koufopanoul (Il994l ). The solid curve shows 



cells, B = (3 /a. Points show data collected by 
the model prediction (1^ for 7 = 1/4 and A = (reproductive cells do not contribute to 
the quantitative trait and the benefit of multicellularity is just a largely arbitrary function 
of the number of somatic cells). An allometric exponent, 7, of 1/5 describes the data best 
(minimizes the mean squared error for data points with i? < 3 x 10^^). The dashed line 
shows the model prediction for 7 = 1/5. The dotted line shows 7^^-B, the limit of gopt for 
5^0. 
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Figure 7: Isolines for organisms of constant size, S = Nr + BNg. The gray curves in the 
lower part of the figure show isohnes with respect to A^^ (see Eq. [2^ . The top part of 
this figure shows the which x corresponds to which A^^.- Isohnes are independent of A and 
approach S/B asymptotically. Decreasing B would increase S/B, change the position of the 
isoline-asymptote and the maximum possible value for x. An increase in A would increase 
^45 and steepen the line that maps a; on A^,.. 
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Figure 8: Optimum number of reproductive cells in organisms with a constant number of 
cells (A^ = 32) that have small somatic cells {B = 0.02, 0.1, 0.2, 0.5, and 1). (a) Isolines (see 
Eq. 1251) change with B (and A) from the purple curve {B = 1) to the gray curve {B = 0). 
I plot isolines for A = 0.1 (solid curves), 0.25, 0.5, and 0.75 (dashed curves). The relation 
between x and Nr is independent of B and shown at the top axis, (b) Fitness of an organism 
as a function of A^^- The curves correspond to the benefit functions from above for A = 0.1. 
For A = 0.25, 0.5, and 0.75 I mark the maximum for each fitness function. Assuming 
that cells are spherical, the figure legend shows the size differences between somatic and 
reproductive cells for the given B values. As expected, decreasing B decreases A^r,opt- The 
position and spread of the maxima shows that for small somatic cells (small B) the ability 
of reproductive cells to contribute to the quantitative trait (parameter A) has little effect on 
A^r.opt and the fitness of the organism. 
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Figure 9: Isolines for organisms of constant size, S = Nr + BNg, that benefit from resources 
that have to be allocated between reproductive cells. The gray curves in the lower part of 
the figure show isolines with respect to A^,. for A = 0.1 (see Eq. 1331) . Isolines change with 
A and B since they contain the term 1/B — A. The top part of this figure shows how A^^ 
chang function of x. This function depends on the parameters A (position of the 

curve) and B (extension of the curve to the right). 
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Symbol Interpretation 

Nr Number of reproductive cells 

Ng Number of somatic cells (an organism is undifferentiated if Ng = 0) 

a Size of a reproductive cell 

P Size of a somatic cell 

B — (5/a The size of a somatic cell relative to the size of a reproductive cell; 

B is bound to [0, 1] 
a Contribution of a reproductive cell to the quantitative trait 

h Contribution of a somatic cell to the quantitative trait 

A — a/h Contribution of a reproductive cell to the quantitative trait relative 

to the contribution of a somatic cell; A is bound to [0, 1] 
X The value of the quantitative trait 

f{x) The benefit function. It captures the extend to which the organism 

benefits from the quantitative trait. The benefit function is scaled 
such that limj;_>oo /(a:) = 1. 

7 The allometric exponent which is approximately 1/4 



Table 1: Notation summary 
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